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Measurement of the Muon {g — 2)-Value 

B. Lee Roberts on behalf of the Muon {g — 2) Collaboration * t 

^'Department of Physics, Boston University, 

590 Commonwealth Avenue, Boston, MA 02215 USA 

The muon {g — 2) experiment is described, and the recent results are presented. These results represent the 
final measurement for the positive muon. 


1. Introduction 

The measurement of magnetic moments has 
been important in advancing our knowledge of 
sub-atomic physics since the famous 1921 paper 
of Stern, Q which laid out the principles of what 
we now call the “Stern-Gerlach experiment”. The 
experimental and theoretical developments in the 
study of the electron’s anomalous magnetic mo¬ 
ment represent one of the great success stories 
of modern physics. The experimental accuracy 
has reached a relative accuracy of ~ 4 parts in 
10® (parts per billion),and the theory is con¬ 
strained by our knowledge of the fine-structure 
constant a, rather than by the eight-order and 
tenth-order QED calculations.^ 


*R.M. Carey, E. Efstathiadis, M.F. Hare, X. Huang, F. 
Krienen, A. Lam, I. Logashenko, J.P. Miller, J. Paley, 
Q. Peng, O. Rind, B.L. Roberts, L.R. Sulak, A. Trofi¬ 
mov - Boston University; G. Bennett, H.N. Brown, G. 
Bunce, G.T. Danby, R. Larsen, Y.Y. Lee, W. Meng, J. 
Mi, W.M. Morse, D. Nikas, C. Ozben, R. Prigl, Y.K. Se- 
mertzidis, D. Warburton - Brookhaven National Labora¬ 
tory; Y. Orlov - Cornell University; K. Jungmann -KVI 
Groningen; A. Grossmann, P. von Walter, G. zu Putlitz 
- University Heidelberg; P.T. Debevec, W. Denninger, F. 
Gray D.W. Hertzog, C. Onderwater, C. Polly M. Sos- 
song, D. Urner, - University of Illinois, Urbana Cham¬ 
paign; A. Yamamoto - KEK B. Bousquet, P. Cushman, 
L. Duong, S. Giron, J. Kindem, 1. Kronkvist, R. McN- 
abb, T. Qian, P. Shagin - University of Minnesota; V.P. 
Druzhinin, G.V. Fedotovich, D. Grigoriev, B.I. Khazin, 
N.M. Ryskulov, Yu.M. Shatunov, E. Solodov - Budker In¬ 
stitute; M. Iwasaki - Tokyo Institute of Technology; H. 
Deng, M. Deile, S.K. Dhawan, F.J.M. Farley M. Grosse- 
Perdekamp, V.W. Hughes, D. Kawall, J. Pretz, S.I. Redin, 
E. Sichtermann, A. Steinmetz - Yale University 
^Supported in part by the USNSF and USDOE 


The gyromagnetic ratio g is defined by 


where s is the spin angular momentum, and g 
is the magnetic moment resulting from this an¬ 
gular momentum. The Dirac equation for point 
particles predicts that g = 2, but radiative cor¬ 
rections increase the value at the part per mil 
level. The magnetic moment g, is defined as 
as g = (1-1- a)efif2m, where a = {g — 2) j2 is 
the anomalous magnetic moment (or simply the 
anomaly). 

When E821 began in the early 1980s, Ou was 
known to 7.3 parts per million (ppm).[Q The 
E821 Collaboration has reported four new mea¬ 
surements with relative accuracies of 13, 5, 1.3 
and 0.7 ppm respectively. 

To the level of the experimental accuracy, the 
electron anomaly can be described by the QED 
of and photons, with the contribution of heav¬ 
ier virtual particles entering at a level of around 
3 ppb. The larger mass of the muon permits 
heavier virtual particles to contribute, and the en¬ 
hancement factor is ^ (m^/we)® 40,000. The 

CERN measurement observed the effect on 
of virtual hadrons at the 8 standard deviation 
level. [Q The standard model value of consists 
of QED, strong interaction and weak radiative 
corrections, and a significant deviation from the 
calculated standard model value would represent 
a signal for non-standard model physics. At this 
conference we have heard much discussion of the 
theory of {g — 2),|9[jl^,|ll|,|2||so it will not be dis¬ 
cussed further in this talk.[|l3[ 
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2. The Experimental Technique 


The method used in the third CERN experi¬ 
ment and the BNL experiment are very similar, 
save the use of direct muon injection into the stor¬ 
age ring which was developed by the E821 collab¬ 
oration. These experiments are based on the fact 
that for g ^2 (or more precisely > 0 ) the spin 
processes faster than the momentum vector when 
a muon travels transversely to a magnetic field. 
The Larmor and Thomas spin-precession and the 
momentum precession frequencies are 
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and the difference frequency gives the frequency 
with which the spin processes relative to the mo¬ 
mentum, 
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which is proportional to the anomaly, rather than 
to the full magnetic moment. A precision mea¬ 
surement of requires precision measurements 
of the precession frequency uja and the magnetic 
field, which is expressed as the free-proton pre¬ 
cession frequency cOp in the storage ring magnetic 
field. 

The muon frequency can be measured as ac¬ 
curately as the counting statistics and detector 
apparatus permit. The design goal for the NMR 
magnetometer and calibration system was a field 
accuracy of about 0.1 ppm. The B which en¬ 
ters in Eq. ^ is the average held seen by the en¬ 
semble of muons in the storage ring, < B > 0 = 
< /o^ Jq^ M (r, 0)B{r, 9)rdrd0 where (j) is 

the azimuthal angle around the ring, r, 9 are the 
coordinates at a single slice of azimuth centered at 
the middle of the 90 mm diameter muon storage 
region. M (r, 9) is the moment (multipole) distri¬ 
bution of the muon beam, and couples multipole 
by multipole with the magnetic held multipoles. 
In the analysis, the held is averaged over the data 
collection time as well. 

The need for vertical focusing implies that a 
gradient held is needed, but the usual magnetic 
gradient used in storage rings is ruled out in our 
case. A sufficient magnetic gradient for vertical 


focusing would spoil the ability to use NMR to 
measure the magnetic held to the necessary ac¬ 
curacy. Furthermore, it is very difficult to obtain 
adequate information on the higher moments of 
the muon distribution in the storage ring, so the 
presence of higher multipoles in < R > is also 
undesirable for this reason. A round beam-prohle 
was chosen, since sharp corners would imply large 
higher moments for M (r, 9). 

An electric quadrupole is used for vertical fo¬ 
cusing, taking advantage of the “magic” 7 = 29.3 
at which an electric held does not contribute to 
the spin motion relative to the momentum. In 
the presence of an electric and a magnetic held, 
the spin diherence frequency is given by 
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which reduces to Eq. ^ in the absence of an elec¬ 
tric held. Note that for muons with 7 = 29.3 in 
an electric held alone, the spin would follow the 
momentum vector. 

The arrangement of a magnetic dipole held 
combined with an electric quadrupole held is 
called a Penning trap in atomic physics. How¬ 
ever with a 14 m diameter and ~ 700 T weight, 
the scale of our trap is quite different from the 
usual one.|H 

In order to meet the conditions discussed 
above, a goal of ±1 ppm uniformity of the 
< B >-held over the storage region was set 
and met. Given the projected knowledge of the 
muon distribution, the allowable strength of the 
quadrupole and higher magnetic multipoles was 
also determined. 

A kick of about 0.1 Tm is needed to bring the 
beam onto a stable orbit. This is achieved with 
three 1.7 m long ferrite-free kickers, fl^ which can 
be thought of as single-loop pulsed magnets car¬ 
rying a current of 4,200 A. The minimum in¬ 
ductance achievable of 1.6 gR limited the peak 
current to 4200A, and resulted in a current-pulse 
width ~ 2.5 times greater than the cyclotron pe¬ 
riod of 150 ns. The phase space mis-match be¬ 
tween the size of the inflector exit and the stor¬ 
age region, and multiple scattering in the inflec¬ 
tor end, reduces the calculated injection efficiency 
to ^ 8.7%. The less than optimal kicker pulse 
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further reduces the injection efficiency to about 
7.3%. Nevertheless the effective data rate per fill 
is almost a factor of 100 over that available in the 
final CERN experiment, Q a factor of 10 coming 
from muon injection, and a factor of 10 coming 
from the AGS intensity. With direct muon injec¬ 
tion, the injection-related background seen by the 
detectors is down by a factor of 50. ||,^ 

The (g — 2) ring functions as a weak focusing 
storage ring with the field index 


n = 


kRq 
PBq ’ 


(5) 


where k is the electric quadrupole gradient. Sev¬ 
eral n - values were used for data acquisition: 
n = 0.137, 0.142 and 0.122, the latter two hav¬ 
ing been used for /i“. The horizontal (radial) and 
vertical betatron frequencies are (approximately) 
given by 


/. = fcVT^ 0.93/c; fy = fcV^ 0.37/c(6) 


where fc is the cyclotron frequency and the nu¬ 
merical values assume n = 0.137. 

The experimental signal is the from fi^ 
decay, which were detected by lead-scintillating 
fiber calorimeters, jr^ The time and energy of 
each event was stored for analysis offline. Muon 
decay is a three-body decay, so the 3.1 GeV 
muons produce a continuum of positrons (elec¬ 
trons) from the end-point energy down. Since the 
highest energy are correlated with the muon 
spin, if one counts high energy as a function 
of time, one gets an exponential from muon de¬ 
cay modulated by the {g — 2) precession. The 
expected form for the positron time spectrum is 

f{t)=Noe~^*[l + Acos{ujJ + (l>)] (7) 


However, a Fourier analysis of the residuals 
from this five parameter fit (see Fig. ^ shows 
a number of frequency components which can 
be understood from the beam dynamics in the 
ring.[^[l8|] The most prominent frequencies in the 
residuals come from the coherent oscillation of the 
beam. 

While the frequencies given in Eq. describe 
the motion of a single beam particle, the aperture 
mis-match between the inflector exit and the stor¬ 
age aperture, along with an imperfect kick and 
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Figure 1. A Fourier Transform of the fit residuals 
from a 5-parameter fit to the 1999 data set. To 
adequately describe the data, additional parame¬ 
ters were necessary (see Eq. H). 


momentum dispersion in the ring, will cause the 
beam to undergo coherent radial oscillations. The 
detector acceptance depends on the radial posi¬ 
tion of the muon when it decays, so any coherent 
radial beam motion will amplitude modulate the 
decay distribution. The principal frequency 
will be the “Coherent Betatron Frequency” 

/cBO = /c - /x = (1 - Vl - n)fc (8) 

which is the frequency a single fixed detector sees 
the beam moving coherently back and forth. This 
motion can be understood by the cartoon in Fig. 
|[ It is this CBO frequency and its sidebands 
from beating with the {g — 2) frequency which are 
clearly visible in the Fourier spectrum of Fig. [|. 
n = 0.137 was chosen to avoid storage ring res¬ 
onances, however the resulting CBO frequency 
was close to the second harmonic of fa{= Wa/27r) 
putting the difference sideband close to fa- The 
n-value was changed for the 2001 data collection 
period to reduce our sensitivity to this difficulty. 

The CBO modified the positron time spectrum 
by 

Np = Noe~{l + A'sin {coat + (j)')) x (9) 

(1 -I- AcBo{t) cos {uJCBOt + (t>CBo)), 
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Figure 2. A cartoon showing the beam motion 
(dashed line), and the apparent motion (solid 
line) seen by a single detector. The vertical axis 
is distance along the storage ring (shown for 4 
turns) and the horizontal axis shows the radial 
betatron amplitude. Ac = “^irp is the cyclotron 
wavelength, is the radial betatron wavelength, 
and p is the central orbit radius. 


3. Results and Conclusions 
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Figure 3. E821 measurements of carried out 
with direct muon injection into the storage ring. 
The relative uncertainties are ±5 ppm (98), ±1.3 
ppm (99), ±0.7 ppm (00), where the number in 
parentheses is the year when the data were col¬ 
lected. All measurements are for /r+. The three 
theory points use the lowest order hadronic con¬ 
tribution from DH98(l^ and the separate DEHZ 
evaluations)^ from r and e+e“ hadronic data. 
The DH98 evaluation includes both r and e+e" 
data. 


where 

A'= A{1 + Ai{t) cos {uicBot + (10) 

4>'= 4>{1 + A 2 {t) cos (ucBot + 4>2)) (11) 

The terms Ai and A 2 can cause artificial shifts in 
uja up to 4 ppm in individual detectors when not 
accounted for in the fitting procedure. The final 
fitting proceedure included these, and other ef¬ 
fects such as pulse pile-up, detector gain changes 
and muon losses from the ring other than by de¬ 
cay. 


After consistent results were obtained in four 
independent (and three different) analysis proce¬ 
dures for uja, and two independent studies of ujp, 
the offsets were removed and the value 


a^+ = 11659 204(7)(5) x 10"^'' (0.7ppm) (12) 


was obtained, [g Excellent agreement was 
found with the previous measurements from 
Brookhaven and from CERN. 

The E821 results are displayed graphically in 
Fig. ^ along with the Davier-H6cker(98)||^ first 
order hadronic contribution, and the most recent 
results from r-decay and e“'"e“ annihilation re¬ 
ported at this conference.110, ll|,[T^,p^ 
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Now that the r-decay value of a^(Had; 1) dis¬ 
agrees with the value obtained from electron 
positron annihilation, the situation is somewhat 
confused. While the path from hadronic electron- 
positron annihilation to a value of a^Had; 1 is 
somewhat more theoretically direct than from 
hadronic tau decay, the extraction of o^(Had; 1) 
from tau-decay data has been carefully studied, 
with all the expected effects included. 

However, as we have heard at this meeting, 
there are systematic differences between other 
quantities when e“'"e“ annihilation and r-decay 
are compared. Since the low-energy e“'"e“ data 
are dominated by the recent precise data from 
Novosibirsk, it is fortunate that further work on 
e"'"e“-annihilation is being carried out at Fras¬ 
cati, Belle and BaBar, and they will either con¬ 
firm or disagree with the recent Novosibirsk data. 
Comparison with the new DEHZ evaluations 
shows either a 1.6 or 3 standard deviation discrep¬ 
ancy. An independent analysis|Q of the e+e” 
data agrees well with the equivalent DEHZ anal¬ 
ysis. Thus the signihcance for an indication of 
new physics will have to wait for clarification of 
the correctness of the hadronic contribution. 

Nevertheless, a recent very conservative evalu¬ 
ation of the impact of {g — 2) on the constraining 
of supersymmetry parameters shows that even 
with the current uncertainties, {g — 2) already 
rules out a “substantial region of (susy) param¬ 
eter space... that has not been probed by any 
previous experiment”. pl|] 

The E821 collaboration has one additional data 
set, which was taken with g,~ in the ring. The 
experimental uncertainty will be on the order of 
0.8 ppm. Obviously we will have to stay tuned for 
further clarihcation on the hadronic contribution, 
and the analysis of the g,~ data. 
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